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Recently, in collaboration with Susskind, we proposed a model of Schwarzschild black 
holes in Matrix theory. A large Schwarzschild black hole is described by a metastable 
bound state of a large number of DO-branes which are held together by a background, 
whose structure has so far been understood only in 8 and 11 dimensions. The Hawking 
radiation proceeds by emission of small clusters of DO-branes. We estimate the Hawking 
rate in the Matrix theory model of Schwarzschild black holes and find agreement with the 
semiclassical rate up to an undetermined numerical coefficient of order 1. 



December 1997 



1. Introduction 



In a recent paper Susskind and the authors presented a model of Schwarzschild 
black holes in Matrix Theory [||. The key feature of the model was the notion of a Boltz- 
mann gas of DO-branes (which we will review briefly below) .0 We showed that the Matrix 
model for M theory in 11 noncompact dimensions, and also for its toroidal compactification 
to 8 noncompact dimensions, contained a set of states which obeyed (up to a numerical 
coefficient which we could not calculate) the Bekenstein-Hawking relation between entropy 
and transverse area for a Schwarzschild black hole. The relation between energy and en- 
tropy for these states also obeyed the black hole formula.@ We were also able to compute 
the long distance Newtonian gravitational force between equal mass black holes, with an 
answer in agreement with classical gravity. 

The purpose of the present note is to calculate the rate of Hawking radiation from our 
model black holes. In we showed that the individual bound DO-branes in the model, had 
the kinematic properties of Hawking radiation in the boosted frame in which we examine 
our black hole. The DO-branes are "tethered" to a classical background by harmonic forces. 
In this note we argue that the probability for the classical variables which produce these 
forces on an individual DO-brane to fluctuate to zero is independent of the mass of the 
black hole in the large mass limit. We show that when this estimate is combined with the 
proper phase space integral, it gives a decay rate for the boosted black hole which is just 
the Lorentz transform of the rest frame Hawking evaporation rate. 

To briefly summarize our black hole model: we consider the Hilbert space of the matrix 
theory which represents Discrete Light Cone Quantized M theory compactified on some 
manifold Y |jl^ in the sector with units of longitudinal momentum {DLCQn)- The 
radius of the light-like circle will be denoted by R. We will also choose the noncompact 
transverse spacetime dimension D to be greater than or equal to 6. The model contains a 
set of variables which includes matrices X* representing the transverse positions of DO- 
branes in a weakly coupled Type IIA string theory called "the analog model." We emphasize 
that these are not the Boltzmann DO-branes of which our black hole is constructed. The 
matrices also describe creation and annihilation operators for strings stretching between 
the DO-branes. 



1 For related work, see §, §, [§ 

^ In the supersymmetric Yang-Mills formulation on the dual torus the corresponding results 
were obtained in 0, g], 
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We consider a semiclassical configuration of the variables of the model, which includes 
a semiclassical background X^^ for the transverse position matrices. This background 
configuration must satisfy a number of constraints, which were described in [|1| . Boltzmann 
DO-brane positions are defined as perturbations of the background 

N 

Xli^Xli + Y,rl,Sr. (1-1) 

n=l 

where the 5n are a set of independent commuting matrices chosen to minimize the quantity 
|Tr[X*;, This is a term in the matrix model Hamiltonian which gives rise to a 

harmonic potential binding these Boltzmann DO-branes to the classical configuration. In 
1^ we argued that for D > 6 we could choose our classical configuration so that this 
harmonic potential did not interfere with the scaling argument we review in the next 
paragraph. 

In matrix theory, the Bose or Fermi statistics of particles arises as a residual gauge 
symmetry. Since the entire gauge symmetry of the model is broken by the classical back- 
ground, this means that the variables should be treated as the coordinates of distin- 
guishable, or Boltzmann, particles. The approximately degenerate configurations of these 
particles then have an entropy of order A^. In we argued that the effective Hamiltonian 
of these degrees of freedom gave rise to bound states whose transverse area is of order GdN 
and whose light cone energy is of order NR{GdN)~'^/^^~'^\ When translated into invari- 
ant mass, this gives precisely the mass/entropy /area relations of a Schwarzschild black 
hole. We were unable to calculate the numerical coefficients in these relations because we 
did not have the full effective Hamiltonian of the Boltzmann gaM and because the mean 
field approximation which we used was very crude. 

According to the above discussion, the light cone energy per particle of the Boltzmann 
gas is ~ R{GdN)~'^/^^~'^\ where Gd is the Newton constant in D non-compact dimen- 
sions. The transverse momentum per particle is ~ [G j:)N)~^/ and the longitudinal 
momentum is just 1/R. We showed that these were precisely the kinematical properties 
of the Hawking particles when boosted into a frame where the black hole has longitudi- 
nal momentum P_ = N/R. This suggests the following attractive and simple picture of 



Recently, Liu and Tseytlin |11] have proposed to describe the interactions of the Boltzmann 



DO-branes by a Hamiltonian containing terms of all orders in the velocity [12|. Perhaps this can 
be used to make some progress on the numerical coefficients. 
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the Hawking evaporation process: the classical background provides a harmonic potential 
which binds the Boltzmann DO-branes to the black hole. The background itself should 
be represented as a coherent quantum state centered around a periodic solution of the 
classical equations of motion of the matrix model. According to the wave function of this 
state, there is a certain amplitude for the part of the classical configuration which interacts 
with a given DO-brane, to fiuctuate to zero. From the point of view of a basis in which 
the particular DO-brane of interest (say di) occupies the lowest right hand corner of the 
matrix, the relevant part of the classical solution is the part which occupies the last row 
and column. There are thus o{N) possible degrees of freedom which might be excited. 

However, in the explicit classical backgrounds which we constructed in D = 11,8 
0], only 0(1) of these possible background degrees of freedom were actually utilized. In 
D = 11, this is a consequence of approximate locality of the action on the world volume of 
the classical membrane. In D = 8 the classical background for a typical black hole could 
be viewed as a lattice of DO-branes connected by strings on a 3-torus, with the lattice 
spacing of order A^~^/"^ Energetic considerations ensured that only strings attaching 
a given DO-brane to o(l) nearest neighbors on this lattice were excited. We will assume 
that similar pictures work for all D > 6. 

As a consequence, it seems reasonable to assume that the amplitude to "liberate a 
DO-brane" from the black hole is independent of as — ^> 00. It is the value of the wave 
function of the system on a submanifold of codimension which is o(l). 



2. Calculation of the Hawking Evaporation Rate 

Given the estimate of the DO-brane liberation amplitude in the previous section, we 
can proceed to calculate the Hawking rate. The quantum fiuctuation of the background 
described in the previous section gives rise to a single DO-brane wave function which is, by 
the estimates of [jl|, a smooth function A{y) of rapid decrease in the variable 

where p± is the transverse momentum and Gd is the D-dimensional Newton constant. 
Thus, the amplitude to produce a DO-brane with momentum much larger than the Hawking 
momentum is highly suppressed. If we assume that the fluctuations which liberate any of 
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the DO-branes are independent and incoherent, then the probabihty per unit time to 
emit a Hawking particle is given by 

where y is the variable defined above. In this equation we have generalized our considera- 
tions to include processes in which a cluster of n DO-branes, with n finite and independent 
of N , is emitted. Such systems also have the kinematic properties of Hawking radiation. 
However, since the cluster of DO-branes is connected to the classical background by o{n) 
degrees of freedom, we should expect the matrix element to fall off exponentially in n. The 
overall factor of in equation (^]^) represents the incoherent sum over processes in which 
a particular bound DO-brane is liberated. Finally, we have used relativistic phase space 
to integrate over the final states of the outgoing DO-brane. Although our calculation is 
done in a particular frame, chosen so that the geometrical structure of the black hole just 
fits inside the DLCQn quantization volume, we expect that the matrix element which we 
have estimated is in fact, for large A^, approximately the invariant S matrix element of a 
Lorentz covariant system. Of course, the crucial, as yet unproven, assumption of Matrix 
Theory is that the S matrix computed by the theory is indeed Lorentz invariant. We 
cannot prove this claim at present, nor can we show that our approximate evaluation of 
the matrix element is accurate. However, assuming that these claims are valid, the correct 
rate is obtained by integrating our matrix element against relativistic phase space. 

Our measure thus contains an extra factor of light cone energy, compared to the 
nonrelativistic phase space of the DO-brane quantum mechanics. This factor, which endows 
the measure with the correct boost transformation properties, has been absorbed into the 
longitudinal momentum delta-function in ( p.2| ) . 

Finally, we need to estimate the square of the matrix element for n DO-branes to be 
liberated, |A(n, y)\'^. As we have explained, this quantity is appreciable only if n and y are 
of order 1. |A(7i,y)p has dimensions of length^"^. The only dimensionful quantities at 
our disposal are R, In and the radii of the compactification torus, Li. Since the measure 
in ( |2.2|) transforms properly under boosts, any dependence of |A(n, y)p on R would violate 
Lorentz invariance. Furthermore, we will assume that the dependence on In and the radii 
is through the Z)-dimensional Newton constant only.S Thus, we are led to 

\A{n,y)\^ r^Go , 

^ This assumption is plausible because Go is the only quantity that appears in the D-brane 
interaction Hamiltonian, but its better justification is clearly necessary. 
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for n and y of order 1. Estimating the Hawking rate (|2.2|) with this assumption, we find 

i?(G,,iV) -2/(^-2) . (2.3) 



By way of comparison, we now compute the Hawking radiation rate according to the 
conventional semiclassical formulae. First, let us write down the Hawking rate in the usual 
equal-time quantization. The answer can be written as 



dN 
d^ 



/ dp+ dp- d^-^p±S{p+p--pl)e-P°/'^''Apo, (2.4) 
Jo Jo J 



where A = 4GdS is the horizon area. We have included the thermal factor appropriate 
for the Boltzmann statistics. The Hawking temperature is related to the Schwarzschild 
radius Rs by 

Th ~ ^ , (2.5) 

and 

Rs ~ (SGd)'^^''-'^ ~ (iVGi>) 1/(^-2) . (2.6) 



An explicit factor of po is needed in (|2.4|) because the measure 



/ dp+ / dp. / rf^-V5(p+P--pi) (2.7) 
Jo Jo J 

is Lorentz invariant. The left-hand side, however, contains a derivative with respect to x'^, 
hence transforms in the same way as po- 

Now we perform a parallel computation in the light-cone frame. Here, the number of 
particles radiated per unit light-cone time is 



dN 
dx+ 



poo poo p 

/ dp- dp+ d^-^p±6{p+p- -p'i)e-P+/'^+e-P~/^~Ap^ 
Jo Jo J 



{2.i 



The factor of is needed for correct boost invariance, since the left-hand side contains a 
derivative with respect to . In the rest frame of the black hole, we have 

yrest _ yrest ^ ~ l/i?5 . (2.9) 

If we carry out a boost 

— R — _|_ Rs -\- 
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then the new temperatures are 



-irest -"^S 1 / D T-irest 



T_ =T'J^^-^ ^ l/R ^ T+ =T1^'^— ^ R/R^s ■ (2-11) 
R Rs 



Doing the integrals, we find 
dN 



^(T+T_)(^-2)/2t+ ~ T+ ~ i?(G'D^) -2/(^-2) . (2.12) 



If we compactify x on a circle of radius R, and work with S ^ N, then the integral 
over p_ is replaced by sum, 



dN 



dx~^ 

n>0 



I V5 - ^) e-^-/^-e-/(^^-^GoN . (2.13) 



While this changes the normalization, the result ( p.l2|) still holds. Thus, our formula ( p.3| ) 
for the rate of emission of DO-branes from one of our matrix theory black holes coincides (up 
to uncalculated numerical factors of order one) with the semiclassical Hawking evaporation 
rate in the light cone frame, ( p. 12 ). 



3. The Rate of Mass Loss 



Using arguments analogous to those that led to ( |2.2| ), we may obtain from the matrix 
model a formula for the rate of light cone energy loss per unit light cone time, 

r dpj d^-'p^5(^-^^)p+\Ain^y)\\ (3.1) 



dx+ 

n>0 

Estimating the integral using the previously stated assumptions, we find 

dE _2 



dx+ 



Now using 



we have 

dM N dE M dN R 



dx+ MRdx+ 2N dx+ i?| 



Boosting back to the rest frame, we obtain 

dM dM dx+ 



l/Rl , (3.5) 



'^^rest d'^^ direst 



(3.2) 



= 2E^ , (3.3) 
R 



where we used the fact that 

dx+ Rs , . 

"•^rest 

The Matrix theory result for the rate of mass loss in the rest frame, ( |3.5D , is consistent 
with the standard semiclassical result. 

What we have shown in this note is that plausible assumptions about the DO-brane 
emission process from the metastable bound state describing the black hole lead to the 
radiation rate consistent with the semiclassical calculations, up to a constant of propor- 
tionality of order 1. Just as in the semiclassical analysis, suppression of the rate for large 
entropy comes from the smallness of the one-particle phase space available at an energy 
comparable to the Hawking temperature [Th scales as A^~^/('^~2)). Our analysis should 
be regarded as a plausibility argument. In particular, we need a microscopic argument 
for why at low momenta the square of the matrix element for liberating a DO-brane is 
of order Go. Nevertheless, we believe that we have described the correct mechanism for 
evaporation of Schwarzschild black holes in Matrix theory. 
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